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Abstract 

We extend our previous work (see arXiv:quant-ph/0501026 ), which compared the predictions 
of quantum electrodynamics concerning radiation reaction with those of the Abraham-Lorentz- 
Dirac theory for a charged particle in linear motion. Specifically, we calculate the predictions 
for the change in position of a charged scalar particle, moving in three-dimensional space, due to 
the effect of radiation reaction in the one-photon-emission process in quantum electrodynamics. 
The scalar particle is assumed to be accelerated for a finite period of time by a three-dimensional 
electromagnetic potential dependent only on one of the spacetime coordinates. We perform this 
calculation in the ^ ^ limit and show that the change in position agrees with that obtained in 
classical electrodynamics with the Lorentz-Dirac force treated as a perturbation. We also show 
for a time-dependent but space-independent electromagnetic potential that the forward-scattering 
amplitude at order does not contribute to the position change in the ^ ^ limit after the mass 
renormalization is taken into account. 

PACS numbers: 03.65.-w, 12.20.-m 
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I. INTRODUCTION 



In classical electrodynamics, when a point charge e with mass m is accelerated by a 4-force 
Fg'^t, it produces radiation and the equation of motion needs to be altered to incorporate 
the radiation-reaction force. In the Abraham-Lorentz-Dirac theory P, 13 we have 

= F!f^, + Fi:^ , (1) 

where are the spacetime coordinates of the charge at proper time r and where the 
Lorentz-Dirac radiation-reaction force is given in the units c = 1 by 
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dr'^ dr \ dr'^ dr"^ 



(2) 



Here we have defined ac = e'^/in, and our metric convention is g^^y = diag (+1, —1, —1, —1). 
We assume in this paper that the external force F^f^^ is a Lorentz force resulting from a 
background electromagnetic field. (See Ref. for an interesting derivation of this equation 
and Ref. ^ for a recent review.) The Abraham-Lorentz-Dirac theory has a number of 
problems if taken without modification, e.g. the existence of unphysical run-away solutions 
in which the charge accelerates under its own radiation. The suppression of this problem 
requires the existence of acausal pre-acceleration. (See, e.g., Refs. 0] for discussion of 
these issues.) However, the treatment of the Lorentz-Dirac force as a perturbation, as 
used and justified in Refs. jsi 1^ for example, is free of these difficulties. Since classical 
electrodynamics is now meaningful only as the h ^ limit of quantum electrodynamics 
(QED), which is the more fundamental theory, it is natural to ask how the perturbation 
theory of QED compares in the classical limit to the Abraham-Lorentz-Dirac theory. Indeed 
we found recently that for a linearly accelerated charged particle the Lorentz-Dirac force 
can be recovered from the one-photon-emission amplitude in QED jiol . 11 , 1^ . In this work 



we studied the position of a charged-scalar wave packet in linear motion and computed 
its change due to the presence of radiation reaction in the h limit. The result was 
found to be in complete agreement with that obtained by using the Lorentz-Dirac force as 
a perturbation in classical electrodynamics. 

In this paper we generalize the work in Refs. [H [ill . Il^ by extending the calculations to 



a charged particle in three-dimensional motion under an external electromagnetic potential 
dependent on one of the spacetime coordinates. We also study the position change due to the 
forward-scattering amplitude at order e^. (Note that Moniz and Sharp derived the Lorentz- 
Dirac force from QED by considering a charge of finite size and taking the zero-size limit . 
Johnson and Hu obtained the Lorentz-Dirac-like force for a classical point charge interacting 



with a massless scalar field by integrating out the scalar field [IJ] . See also Ref. [15| on the 
relation between QED and the Lorentz-Dirac force.) The structure of this paper is as follows. 
Sec. in] outlines the model in the context of classical electrodynamics. We then turn to the 
corresponding model for QED in Sec. Illll and give the expressions for the position shift of 
radiating particle relative to non-radiating particles in terms of the one-photon-emission 
amplitude in the h ^ limit. In Sec.|TV|we calculate the emission amplitudes for potentials 
dependent on one spacetime coordinate and then proceed to calculate the quantum position 
shift in Sec. |3 Then, Sec. IVII gives an expression for the classical position shift due to the 
addition of a radiation-reaction force treated as a perturbation. This quantity is found to be 
identical with the h ^ limit of the quantum position shift. We summarize and discuss our 
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results in Sec. IVIII In Appendix |X] we demonstrate that the forward-scattering amphtude 
at order does not affect the position of the particle in the h ^ limit after the mass 
renormalization is performed. 

II. EXPLANATION OF THE MODEL 

In this section we shall describe the model used for the investigation in this article in the 
context of classical electrodynamics. The model is the three-dimensional extension of that 
used in Refs. 3| with the potential dependent only on one of the space-time coordinates, 



x"" say. The particle is accelerated by an external force arising from an electromagnetic 
potential V^. That is, the external force in Eq. is given by 

n T 

F!:^, = {d^V'^-dWn-^. (3) 

We wish to calculate the position shift, i.e. the change in the position of the accelerated 
charge due to radiation reaction, after a finite period of acceleration. Thus, we shall assume 
that the potential behaves as follows: V^{x'^) = V^^^, where V^g-^ is a constant vector, 
for x°' < —Xi and V^{x°') = for x°' > —X2 with Xi and X2 being positive constants 
{Xi > X2). Hence, the acceleration occurs only in the region — Xi < x"' < —X2. We also 
require that the coordinate of the particle increases as a function of time. (This condition 
is automatically satisfied if x° is the time coordinate.) For the sake of simplicity we set up the 
model so that the charged particle would pass the space origin at time t = in the absence 
of radiation reaction. (The spacetime origin is singled out as a convenient point of reference 
because the quantum fields are expanded in terms of the functions e^*P'^/^, which take the 
value one at the spacetime origin. See Ref. for an analysis of the model in one dimension 
with a reference point other than the origin.) The conditions on the external electromagnetic 
potential imply that at t = the particle has already gone through the region of acceleration. 
The position shift is defined to be the change in the space coordinates {i = 1, 2, 3) of this 
particle away from = at t = due to radiation reaction. 

It will be convenient for later use to write down the space components of the Lorentz- 
Dirac force explicitly in terms of the velocity = and acceleration a* = u*, where the 
dot represents differentiation with respect to the coordinate time t. We write the equation 
of motion for the Abraham-Lorentz-Dirac theory ([Q) in a slightly different form as 

d ( dx^\ , ,,,dT ■ 

Then, defining x, v and a to be the vectors with i-th components x*, and a* respec- 
tively, and introducing the definition 7 = dtjdr = (1 — v^)~^/^ as usual, we obtain after a 
straightforward calculation 

= {I [7'(^ ■ + 7^«1 - 7^(^ ■ V) V - . (5) 

III. POSITION-SHIFT FORMULA IN QED 

In this section we present the expression for the h —>■ limit of the quantum position shift 
in terms of the one-photon-emission amplitude. The Lagrangian density for the quantum- 
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field-theoretic model of a charged scalar field with mass m and charge e coupled to the 
electromagnetic field in the Feynman gauge is given by 

C = [{D^ + ieAJh)^]\D^^ + leA^/n)^ - {m/Uf^^^ - ^F^.F^'^ - \{d^A^f , (6) 

where F^^ = d^A^ — dyA^ and = + iV^/h. The function is the external potential 
which accelerates the charged scalar particle as in the previous section. For e = the 
charged scalar field is expanded as 

^^"^^ = ^ / 2po(2^^)^ [^(P)*p(^) + B\p)%{x)] (7) 
with po = y/p^ + m^, where the non-zero commutators are 

[A{p), At(p')] = [Bip), i?t(p')] = 2po(27r;.)353(p _ p') . (g) 

The mode functions $p(a;) satisfy the Klein-Gordon equation in the presence of the external 
field V^: 

(Z}^D^ + m2)$p(a;) = 0. (9) 

If the potential is purely t-dependent, then we require that $p(a;) take the form e"*^'^/^, 
where p-x = p^t — p ■ x, near the hypersurface t = 0. (This is possible because the potential 
Vf^ vanishes there by assumption.) If the potential is space-dependent, then we require that 
$p(x) take the form ^-^v-^lft. a region including the spacetime origin Similarly, the 
electromagnetic field is expanded as 

^^^K(k)e-*- + at(k)e*^] (10) 

with k = ||k|| and A; ■ x = /ct — k ■ x, where the annihilation and creation operators, a/i(k) 
and aji(k), for the photons with momenta hh satisfy 

[a,{k), 4 (k')] = -g^,i2nf2hk5%k - k') . (11) 

Notice that the scalar field ip is expanded in terms of the momentum p whereas the elec- 
tromagnetic field Afj^ is expanded in terms of the wave number k. We adopt this convention 
because the vectors p and k are regarded as classical rather than p/h, the wave number of 
the scalar particle, and ^k, the momentum of the electromagnetic field. 
Let the initial wave-packet state be given by 

where the function /(p) is sharply peaked about a given momentum with width of order h. 
If the potential depends on a space coordinate x°', then we require that the a-th component 
of this momentum be positive so that the wave packet is near the spacetime origin after 
it has undergone acceleration. The normalization of the operators A^(p) is such that the 
condition = 1 leads to 
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(12) 



This shows that the function /(p) can heuristically be regarded as the one-particle wave 
function in the momentum representation. The position expectation value is given by 



{x^)= j rf3xx^(p(x,t)), (14) 

where the density operator for the scalar field is 

p{x) = ^-.ip^'dtip - dt(p^ ■ (p: . (15) 

In Ref. it was shown that for the non-radiating particle the position expectation value 
at i = is given to lowest order by the expression 



where dpi = dpi — dp'- This expression can heuristically be regarded as the expectation 
value of the position operator ihdpi in the momentum representation. It also agrees with 
the expectation value of the Newton- Wigner operator 3]. We let the function /(p) satisfy 



(a;*)o = so that we measure the position shift for the particle relative to the origin as we 
required in the previous section. 

For the radiating particle, the evolution to order in perturbation theory in interaction 
picture takes the following form: 

At(p)|0)-[l+.^(p)]At(p)|0) + - j ^^^^^(p,k)at(k)At(P)|0), (17) 

where J^(p) is the forward-scattering amplitude and where w4'^(p, k) is the amplitude for the 
one-photon emission. If the potential is purely t-dependent, then the momentum is conserved 
and one has P = p — /ik. If the potential is x^-dependent, say, then the momentum P is 
determined by energy conservation -^p^ + m? = y/P^~+^v? + hk and transverse- momentum 
conservation = P* + hk^ [i = 1,2). Thus, the final state for the initial wave-packet state 
1^) given by Eq. (|T2|l is 

The position expectation value for this state in the ^ ^ limit can be calculated by using 
Eq. (fT^ and was given in Ref. jl^l- It consists of three terms identified as the non-radiating 
position expectation value, (x')o, which we have set to zero for simplicity, the forward- 
scattering contribution (which is not due to the radiation) and a third term identified as 
the position shift due to radiation reaction, which we denote by Thus, the position 

expectation value at t = is 



(x^)<=o = -/^^Re^(p) + 5a;^ (19) 
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with 

Sxh = -2 y ^^(^^'*(P' k) AiP, k) (20) 

in the ^ — ^> hmit, where the wave packet becomes narrowly peaked about the momentum 
p. (From now on, the momentum p will always refer to the peak value for the wave-packet 
state.) Thus, the position shift 6xq can be found from the emission amplitude ^^(p,k) 
to be calculated for a given external electromagnetic potential that accelerates the particle. 
Although this expression was considered only for linear acceleration in Ref. [12], it is equally 
valid for motion in three dimensions: the details of the potential and the path of the particle 
are unnecessary for its derivation. 

It will be shown in Appendix A that the forward-scattering contribution in Eq. (|19|) 
vanishes in the ^ — >■ limit after the mass renormalization is taken into account if the 
potential depends only on time, and we expect that this is also the case for potentials 
dependent on one space coordinate. For this reason we concentrate on the evaluation of the 
radiation-reaction contribution 6xq in the rest of this paper. 



IV. EMISSION AMPLITUDE 

The photon-emission part of the evolution of the state in Eq. (fT7|) is given in terms of 
the emission amplitude by 

^t(p)|0)^... + -y ^^^^^(p,k)at(k)At(P)|0). (21) 

The evolution of this state to first order in time-dependent perturbation theory is 

^t(p)|o) - ■ ■ ■ - ^ 1 d'xnjix)A\p)\0) , (22) 

where 'Hi{x) is the interaction Hamiltonian density. Comparing these two evolution expres- 
sions, we can write the emission amplitude in terms of the interaction Hamiltonian density 

as 

A(P,k) = \j j d'x{Q\a,{}.)A{jp')'Hr{x)A\Y>m . (23) 

The Hamiltonian density can be obtained from the Lagrangian density © by the standard 
procedure. Thus, we find the following interaction Hamiltonian density: 

Hj{x) = -A,:[¥,t/)'^v'- (Z}^V)V]: + ^5^/l,A,:y,ty,:-_^:y,V. (24) 

i=l 

where = d^ + iVf^/h as before. We have normal-ordered the scalar- field operators to drop 
the vacuum polarization diagram automatically. (Note that the second term is different 
from what might be naively expected, —(e'^/fi^)A^A^:(p'^ip:. This difference is due to the 
presence of interaction terms involving (f or (f^^ in the Lagrangian density.) The last term in 
Eq. fj24|) is the mass counterterm. The counterterm for the wave-function renormalization 
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will not be necessary in our calculations. A photon is emitted due to the first term in the 
Hamiltonian density ()24|). Thus, we have 

^,(p,k) = ^ I ^^j^, J rf^x(0|a,(k)A(p'){A.:[v,t/)V-(^V)V]:}^nP)|0). 

(25) 

By using the expansion of the fields and ip, and the commutation relations for the 
annihilation and creation operators, one readily finds 

^,(p,k) = -^eh I I d'xe''- {%,{x)D,%ix) - [D,%,ix)]^ %ix)} .(26) 

Next we shall find a simple expression for the h ^ limit of the emission amplitude 
given by Eq. for an electromagnetic potential which depends only on one spacetime 
coordinate. We present the details of the calculation only for a t-dependent potential. The 
final result for an x- y- or z-dependent potential is identical, and the techniques involved 
are mostly the same. Thus, we consider a time-dependent potential V'^{t) {i = 1,2,3) with 
V^{t) = 0. [If V^°(i) 7^ 0, one can gauge away this component.] The system is translationally 
invariant in the spatial directions and, hence, we can let 

$p (x, t ) = 0p (t ) exp (zp ■ x/;i) . (27) 

The amplitude in a spatial direction for the t-dependent potential is then 

A\p,k) = -ej I rf^x0;,(t)0p(t)[p^ + /-2r(t)]e^[(P~P')-/^le^('=*-'^-) 

= -e [ rfte'^Vp(^)0p(^)^-^^^ (28) 
J Po 

with P = p — hk, where we have let P* +p* — 2V'^{t) = 2[p* — V{t)] as the difference — 
is of order h jl^]. It would be wrong to equate 0p(t) with 4>p(t) because these functions 
oscillate with periods of order h~^. For the time component we have Dt = dt and simply 
obtain 

A'ip, k) = / dt [05^5*0p - {dt<P*p) 0p] e''' . (29) 

To proceed further we need to approximate the function (f)p{t) in a way suitable for taking 
the h ^ limit. To this end we use the semi-classical WKB approximation. By substituting 
Eq. (P7|) in the Klein-Gordon equation Q we find 

{h^d^ + [p - Y{t)f + m^} 0p(t) = , (30) 

where V(t) is the vector with its z-th component given by V^{t). The standard WKB 
approximation gives the following positive-frequency solution: 



, ^p(^) 

where 



(31) 



ap{t)^J[p-Y{t)Y + m^. (32) 



We note that the local momentum and energy of the point particle corresponding to the 
wave packet considered here are 

= P - , (33) 
= a^it). (34) 
Now, the product of two wave functions in the emission amplitude ()28|) can be written 

<pmMt) = ^e^p {-J- 1 [^p(o - ^p(c)] ^c} , (35) 

where we have replaced Pq and crp(t) in the pre- factor by po and crp(t), respectively, because 
we are interested only in the h ^ limit. The integrand in the exponent can be evaluated 
to lowest order in h by using Eqs. (jH^ and (jH^ as 

- (36) 

where the repeated indices i are summed over. By substituting this approximation in 
Eq. dnS)) we find 



= — exp (— zk • x) , (37) 

where we have used the fact that the particle passes through the spacetime origin. By 
substituting this formula in Eq. (PH|1 and noting Eqs. and we obtain 



d^^e''^, (38) 

where we have defined ^ = t — n • x with n = k/Zc. We emphasize that x and ^ here 
are functions of t evaluated on the world line of the corresponding classical particle passing 
through the spacetime origin. 

Let us now consider the time component ^°(p, k) of the emission amplitude given by 
Eq. (j2ni)- Note that from the WKB expression for 0p(t) we have, to lowest order in h, 

dtMt) = -^^p(^)0pW • (39) 
By substituting this formula in Eq. we obtain 

A%p,k) = [dt Mt) + Mt)]<put)Mty 

'^Po J 
-e / c/te-^'^-V'^* 

^^e^'' , (40) 
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where we have let ap(t) = ap(t) and used Eq. (|H7j) . By combining this formula and Eq. (|HH)) 
we obtain the following concise expression for the ^ ^ limit of the emission amplitude: 

^e^'^ . (41) 

One can proceed in a similar manner for the case where the potential depends on any one of 
the spatial coordinates. One obtains Eq. ()4ip in this case as well. The emission amplitude 
fl41|) is identical to that for a classical particle corresponding to the wave packet considered 
here 



V. QUANTUM POSITION SHIFT 

The expression ()4H1 for the emission amplitude is ill-defined, as noted in Ref. for 
the one-dimensional case, because the integrand does not tend to zero as ^ — *• ±oo. For 
this reason we introduce a smooth cut-off function xiO which takes the value one while the 
acceleration is non-zero and has the property limg_^-|-oo x(0 = 0, viz. 

A^{p,k) = -eJ d^^xiOe''^. (42) 

One can show exactly as in Ref. 112] that the substitution of this expression into the quantum 
position shift given by Eq. ^2U^ leads to an expression independent of the cut-off function 
as follows: 

d'^x,,, d ( dx^ \ 



where dfl is the solid angle in the wave-number space of the photon emitted. The variable 
appearing as the subscript for a partial derivative — the variable ^ in this equation — is 
held fixed. [Note that dxf'/dt = {dx^'/dt)^,] 

To evaluate (5xq in Eq. (pSj) we first need to put the integrand in the form amenable to 
the solid-angle integration by rewriting it in terms of t instead of ^. One can readily write 
d'^x'^/dC,'^ in terms of t-derivatives by using d/dC, = (1 — n'^x'^)~^d/dt as follows: 

^ = ^-3 r (1 _ n'x') x^ + n'-£x^\ , (44) 

where ^ = 1 — Here and in the rest of this section, Latin indices take the spatial values 
1 to 3, and are summed over when repeated. The time and space components of Eq. ()44|1 
can separately be given as 



di^ 

d^x^ 



^=r'n'x\ (45) 
[(1 - n'x') F + n'x'x^] . (46) 



Next we express (dx^/dp^)^ in Eq. ()43p in the form involving t rather than ^ as follows. Note 
first 
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[The zeroth component of this equation is trivial because {dt/dp^)t = 0.] By substituting 
dt = + n^dx^ in this equation with fi = j and solving for rfx*, we obtain 



where = x^. Hence 



V 



1 — n^v^ 



d^ + 



[5*^(1 — n^v'') + n'^v'^] f dx 



1 — n'-v 



lr)d 



dpi 



dp' 



dx^ 
dp^ 



[6 {1 — n V ) + n f *] / dx 



1 — n'f ' 



dp^ 



(48) 



(49) 



With ^ fixed we have dt — n^dx"^ = dC, = 0. Thus, 



dph> J 



dp^ 



1 — Tn}v^ ydpi J t 



(50) 



By substituting Eqs. (jlKjl . (jlEI), P^j) and (jHUjl in Eq. we find the following expression 
after a tedious but straightforward calculation: 



ll'T^a^ - ha^ - lia'v^ - /^^(a ■ v) 
if-i-^a^a' - 2I^'a^{si ■ v) - jfa^ 



1- 


f dx''\ 


J rft 


^df ) 



dx'' 
dp^ 



where 



-'2 



jijk 
-'3 



^2 ' 

^4 3 ' 3 ' 



dQ- 



(51) 

(52) 
(53) 
(54) 
(55) 



Evaluation of these solid-angle integrals is facilitated by noting that the last three integrals 
are proportional to partial derivatives of Jq with respect to v\ Substitution of Eqs. (j3^ -(fS ^ 
in Eq. yields 



6x 



2ar 



Q 



(56) 

By comparing this equation with the expression (jS)) for the Lorentz-Dirac force, we note 
that the quantum position shift can be written as 



Sxl 



dx^ 
dp^ 



(57) 
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We have used the fact that a(t) 7^ only for t < and for a finite interval of time to 
integrate the first term in Eq. (I^Uj) by parts. In the next section we shall demonstrate that 
this expression for the position shift agrees with that in the Abraham-Lorentz-Dirac theory. 
We recall here that is the momentum of the particle at t = 0, which is after it has gone 
through the acceleration. 



VI. CLASSICAL POSITION SHIFT 

In this section we discuss the classical perturbative solution to the Lorentz-Dirac equation 
(0)) in general terms in order to show that the classical and quantum position shifts coincide. 
We do not need the assumption that the external electromagnetic potential depends only 
on one coordinate in this section; this assumption was necessary only for the calculation of 
the quantum position shift Sxq. 

Recall that the classical motion of a charge e with mass m in an electromagnetic potential 
= (y^, V), which depends on t and x, is described by the following Hamiltonian: 



H = ^(P -Vy + m'^ + eV\ (58) 

where P is the momentum conjugate to x. One can readily show that the Lorentz-Dirac 
equation (jD) is equivalent to the following set of equations: 

f) J-f 

^^ = -^+-^LD- (60) 

Now, let (x, P) = (xo(t), Po(t)) be the solution to these coupled equations with jF-^j-, = 
satisfying xo(0) = and Po(0) = p. This solution gives the classical trajectory of the 
particle passing through the spacetime origin with final momentum p in the absence of 
radiation reaction. If (x, P) = (xo(t) +5x(t), Po{t) +6P{t)) is the retarded solution to these 
equations with jF^j-, 7^ to first order in J-'l^, then 5x and 6P will have the property that 
((5x, 6P) (0, 0) as t ^ —00 and satisfy 

J pp. TT f)2 TT 

where the partial derivatives of H are evaluated at (x, P) = (xo(t), Po(t)) and, hence, are 
functions of t alone. The quantity 6xq = 6x^(0) is identified as the position shift in the 
classical Abraham-Lorentz-Dirac theory. Thus, our task is to show that 6xq = Sxq. To this 
end we define a set of solutions to the homogeneous equations obtained by letting JF^^, = in 
Eqs. dnH) and Thus, we define the solutions {6x'{t),6P\t)) = (AxJ^.)(t; s) , AP'^.^{t; s)), 
with labels j = 1, 2, 3 and s G (—00, +00), to these homogeneous equations by the following 
initial conditions: 

Axys;s) = 0, (63) 
AP^)(.;s) = 6'^. (64) 
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Then, the solution (Ax*^^-j(t; s), AP^*^.-|(t; s)) represents the particle trajectory which coincides 
with xo(t) at t = s and has an excess momentum solely in the j-direction at this value of t. 
Then it can readily be seen that the retarded solution to the inhomogeneous equations 
and (jH^ is given by 

6x' = I dsJ^lj^{s)Ax\^-^{t-s) , (65) 

J — oo 

6P' = f dsJ'Us)APi^^{t;s), (66) 

J — oo 

where the index j is summed over. Hence, the classical position shift is 

6xh= f dtTUt)Ax\^^{Q-t). (67) 

J — oo 

To obtain a similar expression for the quantum position shift 5x^q we note first that 

' dx^ 



dp' 



A4)(t;0). (68) 



This is because eAx'|.^(t; 0) with e -C 1 is the change in the position of the charged particle 
in the j-direction at time t caused by a change in the i-th component of the momentum at 
t = by e. Thus, the quantum position shift Sxq in Eq. can be written 



6x'q = -[ dtJ^Ut)Axl^{t-0), (69) 

J — oo 

where the index j is summed over. 

The classical and quantum position shifts given by Eqs. (jHTj) and (j69|l . respectively, can 
be shown to be equal by using conservation of the symplectic product for the homogeneous 
equations (jHT|) and (jH^ with jF-^p = 0. The symplectic product between two solutions, 
(Ax^, AP^) and (Ax^, AP^) to these homo geneous equations is defined by Ax^ ■ AP'^ — 
Ax'^ ■ AP^ . It is well known and can easily be verified that the time derivative of this 
product vanishes, i.e. that it is conserved. By equating the symplectic products of the two 
solutions (Aa;^.)(t; s), AP^^-^{t] s)) and {Ax'^j^{t; u), AP^.-^{t] u)) aX t = s and u we have 

Ax\^ {s- .) AP(^.) {s- «) - A4) {s- u)AP^^ {s-s) = Ax\^ {u- s) A^) {u- u)-Ax\^^ {u- u)AP^^ {u;s). 

(70) 

This equation and the initial conditions and (jMj) imply 

-AxJ^.)(s;m)=AxJ,)(m;s). (71) 

In particular, we have —Axl-JO; t) = Axj.Jt; 0). This equality and Eqs. (jU7j) and (jUU]) imply 



that 5a;^ = 6x1 
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VII. CONCLUSION 



In this paper we compared the change in position of a charged scalar particle due to radi- 
ation reaction in both classical and quantum electrodynamics in three space dimensions. We 
found that, for a charged particle accelerated for a finite period of time by an electromagnetic 
potential dependent on one of the space time coordinates, the change in position in classical 
electrodynamics using the Lorentz-Dirac force as a perturbation is the same as that given by 
the h ^ limit of the one-photon-emission calculation in quantum electrodynamics. This 
extended the results of Ref. 2| from one to three dimensions. Since the forward-scattering 
amplitude at order does not contribute to the change in the position in the ^ ^ limit 
if the potential depends only on time as shown in Appendix A, one can conclude that the 
^ ^ limit of quantum field theory coincides with the Abraham-Lorentz-Dirac theory at 
least in this case. We expect that this will be true also with a potential dependent on one 
space coordinate. 

In our calculations it was crucial that a plane wave corresponded to a set of parallel 
classical trajectories with the same motion. This property was guaranteed by the assumption 
that the external electromagnetic potential depended only on one spacetime coordinate. If 
the potential depends on more than one coordinate, classical particles with the same initial 
or final momentum can have many different motions depending on their initial or final 
positions. (Consider, for example, scattering by a central potential.) It would be interesting 
if one could overcome this difficulty and extend our results to cases with more general 
potentials. It would also be interesting to extend our results to radiation reaction in curved 
spacetime or to the back-reaction to gravitational radiation 



APPENDIX A: FORWARD-SCATTERING AMPLITUDE 

In this Appendix we show that the contribution of the forward-scattering amplitude at 
order to the position change of the particle vanishes in the ^ — > limit if the external 
electromagnetic potential depends only on t. Since this contribution has an explicit factor 
of h in Eq. ()19p. we can neglect any term of order higher than in the forward-scattering 
amplitude. What we demonstrate here is that the one-loop contribution of order or 
lower is cancelled exactly by that from the mass counterterm, which is formally of order 

We use a non-covariant Hamiltonian formulation with the Hamiltonian density (j24|) . The 
forward-scattering amplitude J-'{p) for the state with final momentum p is given by 

L I d^x'd^x{Q\A{p')T[ni{x')ni{x)]A\p)\Q) , (Al) 

where T denotes time-ordering, to second order in the standard time-dependent perturbation 
theory in the interaction picture. 

We first consider the contribution of the mass counterterm —[5rn?/Ti'^) : ip'^ip : in the 
interaction Hamiltonian density to the forward-scattering amplitude. This contribution. 
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which we denote by jF™'^***^(p), can readily be found at first order in perturbation theory as 

. . 1 



2%o 



(A2) 



where 0p(t) is the time-dependent part of the scalar mode function defined by Eq. (j2H). 
(The quantity jF™'^'^'^(p) is obviously divergent in the t-integration. In the end all terms with 
this divergence property will cancel out.) The mass parameter in the counterterm in the 
standard covariant perturbation theory reads 



= — 

h 



{27c)H { [g2 - m2 + ie] [{p - qf + it] [{p - qf + ie] 



{p + qf 



(A3) 



It is convenient to integrate over the time component go of for later purposes. The result 
can be given as follows: 



5m 



q 



n J (27r)^ 



(p + q) 

^Kqo 
3 1 

+ 77T7 + 



1 



2K 4Kqo 



qo + K-po K + po + qo 
(po-go)^ ^ (Po + qo) 



qo + K + po qo + K-po 



(A4) 



with K = p — q and K = ||K||, where go is now defined to be \/c[^ + 



The second term in Eq. ()24|) contributes to the forward-scattering amplitude at first order 
in time-dependent perturbation theory. By denoting this contribution by J^2{p) we have 



i=l 



2tpo{2nnfj^2{p)S'{p - P) = E / {0\A{p)A{x)A{xy.v^{xMxy.A^{p')\0) . 

(A5) 



One can readily evaluate this expression with the following result: 

3e2 f+°^ , f d^q 



2h^Po 



2K{2-Kf ' 



(A6) 



where we have used d^\i/k = h~'^d'^q/K. 

The first term in Eq. contributes to the forward-scattering amplitude in second order 
in time-dependent perturbation theory. Writing this contribution as J^i(p), we find 



-^i(p) 



d' 



q 



2j9o J 2qo{2TTn)^ 2K 



dtidt2 

e{tl-t2) [0;(tl)0p(t2) Bl(tl,t2,P,q)0q(tl)0;(t2) 
+ e{t2 - h) U;(ti)0p(t2) Vl {tl,t2, P, q)0;(tl)0q(t2 



-iK{h-t2)/h 



jK{ti-t2)/h\ 



where 



Pi(ti,t2,P,q) = dt.dt, + [p + q - 2V(ti)] • [p + q-2V(t2)] • 



(A7) 



(A8) 



It is convenient to define the variables t and rj hj ti = t — h'q/2 and t2 = t + hT]/2. Then 



•^i(p) 



2h^Po 



2go(27r)3 2K 
14 



rft[G'_(p,q,t) + G'+(p,q,t)] 



(A9) 



with 



G_(p,q,t) 
G+(p,q,t) 



dr] [0;(ti)0p(t2) 252(t,r/,p,q)0q(ti)0;(t2; 

dr] U;(ti)0p(t2) 52(t,^,P,q)0;(tl)0q(t2) 



-iKri 



where 



B2 (t, r/, p, q) = -- - + [p + q - 2V(ti)] ■ [p + q - 2V(t2) 



(AlO) 
(All) 

(A12) 



We write the time-dependent part of the mode function, 0q(t), as 



(A13) 



where '?/'q(t) — > 1 as /i — > 0. By substituting this expression in Eq. (jHUIl we can expand ipc^if) 
in terms of Ti as 

,lj^{t) = l + ihij'^^\t) + 0{h^). (A14) 



The exphcit form of T/^q ^ (t) is unnecessary though it can easily be found. Note also that 

r.7y/2 
-r,/2 



0q(^l)0q(^2) = ^====exp 



r-ri/2 



It can readily be shown that the functions ^^(p, q, t) are of the form 



(A15) 



/•±oo 

G±(p, q, t) = ± / rfr/ [/±(p, q, t) + ^(/i^)] 



X exp 



where the function /±(p, q, t) can be found as 

/±(P, q, = {- WAt) T cr^{t)f + [p + q - 2V{t)f} |0p(t) H^qW P • 



(A16) 



(A17) 



An important fact to note here is that there are no terms of order h in the pre-factor inside 
the first square brackets in Eq. ()A16|) . This fact is a simple consequence of the equations 
rciti)Mh) = 1 + 0{h^) and 



[p + q - V(ti)] ■ [p + q - V{h)] = [p + q - V{t)f + 0{h') 



(A18) 



Let us first consider the integral ^^(p, q, t). We change the integration variable from rj 
to (3 defined by the following relation: 



[Orp(t) + Orq(t) + K](3 = / [dp (t + kQ + (Tq(t + kQ + K] dC . 

J-n/2 



(A19) 
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Expanding the integrand and integrating, we solve for as a function of (3 for small h and 
find 

1 ap(t)+<Tq(t) _f^2p2^(j^^ApA^ 



and 



Then we have 



r/ 



dri 



1 



24(Tp(t)+(Tq(t)+ir 



(A20) 
(A21) 



/■OO 

G+(p,q,t)= / rf/3[/+(p,q,t) + 0(^2)] exp{-^[ap(t) + aq(t) + ir]/3} . (A22) 
Jo 



Introducing a convergence factor by replacing K hy K — ie, we obtain 

«/+(p,q, 



G+(p,q,t) 



aq(t) + ap(t) + 

The corresponding term in the forward-scattering amplitude is 

d^<i 1 /+(p,q,i) 



(A23) 



•^i+(p) 



2^'po 



(it 



2go(27r)3 2ir (jp{t) + a^{t) + K 



(A24) 



The amplitude jFi^(p) can readily be seen to be ultra-violet divergent. 

Next we analyze the contribution from G^{p, q, t). One cannot proceed as above because 
of the infrared divergence in the q-integration as we shall see. We define the variable (3 in 
analogy with the variable f3 in Eq. ()A19|) by 



r,/2 



ri/2 



-ap{t + hC) + aq(t + hC) + K]dC. 



Now, for small K = ||p — q||, we have 

-fTp(t) + (Tq(t) +K^K- v(t) • K , 



(A25) 



(A26) 



where v(t) = [p — V(t)]/o"p(t) is the velocity of the classical particle with final momentum 
p (see Eq. (I36|) ). Hence, in the limit K one finds 



1 



1 - v(t) ■ n 



v/2 



v/2 



[l-^(t + hC)-n]dC., 



(A27) 



where n = K/K. Thus, if we write dr] = J{p,q^,t,hP)dp, then the function J(p,q, 

is finite as — ^ 0. Hence, the expression corresponding to Eq. (jA22j) can be given in the 

following form: 



G_(p,q,t) 



d/3 



/_(p,q,t) + 5^;i"/3^/„,_(p,q,t) 



n,d 



-«/-(p,q,^) 



i)''(i!r/,,-(p,q.t) 

d+1 



ap(t) + aq(t) + K ^ [-ap(t) + aq(t) + K] 



(A28) 
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with n > 2 and n > d > 0, where /„d_(p, q, t) are finite as 0. 

Now, if we substitute Eq. ()A28|) in Eq. ()A9|) . then the terms with d > 1 are infrared 
divergent in the q-integration because \imK^o[—o'p(t) + (Jq{t) + K] — >• as can be seen from 
Eq. ()A26|) . (Note that d^q = d^K..) To deal with this difficulty we cut off the integral over 
q by requiring K > Kq = h"\ with | < a < 1, where A is a positive constant, and postpone 
the analysis of the contribution from K < Kq till later. (The reasoning for this particular 
choice of limits for a will be seen later). Then, we find that the small- i^' contribution of 
each term in Eq. (fX28|l to the q-integral behaves like H'^Kq-'^ = ff+^^~'^)'^\^-'^ if d > 2 and 

log(^"A) ii d = 1. Now, since 1 — q;>0, n>2 and n > d, we have n + (1 — d)a >2 — a. 
Therefore, in the h ^ limit, the fnd- terms will not contribute above the cut-off, leaving 
only the first term. Let us thus combine the leading order terms from J-'i (using and the 
first term of G_, above and below the cut-off) and the result for JF2 to define 

/i-V f , , ,:2 f d\ 1 



2crq(t)(27r)3 2ir 



, -[ap(t)-aq(t)f + [p + q-2V(t)f ; 
+ a^{t) + K + a^{t) 

(For f± we have used the relation |0qp = qo/(J^{t) \l + 0{h'^)\). What remains is the 
contribution of minus the first term (which is in jF'^^'^^s)^ from below the cut-off. Thus 
we define 

2po JK<h'-x 2qo{27rhy 2K J 

X0;(ti)0p(t2) 5i(ti,t2,P,q)0q(ti)0*(t2) (A30) 



•^<'°(P) = ^ [dt\Mt)\' [ 
2po J Jr 



K<n^^ 2aq(t)(27r)3 2ir 
^^ -[ap(t) + aq(t)f + [p + q-2V(t)f 

a^it) + K - ap{t) ' ^ ^ 

with the operator "Di (ti, ^2, P, q) defined by Eq. ()A8j) . Hence the forward-scattering ampli- 
tude can now be written as follows: 

^(p) = ^--«(p) + ^i<=^<ii'^g(p) + ^<(p) - ^<'°(p) + 0{h-") , (A32) 

with the contribution from the mass counterterm, jF™'*^^(p), given by Eq. ()A2|) . Below we 
show the equality jF'°'**^™^(p) = — jF™'^'^'^(p) and demonstrate that the quantity jF<(p) — 
jF<'*'(p) is of order and is purely imaginary at this order. Then one can conclude that 
the real part of the forward-scattering amplitude is of order higher than h~^, and hence does 
not contribute to the position change, as claimed in this paper. 

The equahty ^^*^'^*^™s(p) = —J^^^^(p) can readily be obtained by letting q = q— V(t) and 
p = p — V(t) in Eq. ()A29|) . Then writing o"p(t) = ^/^~+~rn? = po and (Jait) = \/^~+~rn? = 
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go, we have 



pleading 



2^2 



V) 



'^0 



2go(27r)3 2fs: 



X 



, -(po + go) +(p + q)^ , -{po-qo 

-6go H . — z h 



(p + q)' 



qo + K -po 



qo + K + po 



(A33) 



This can be seen to be the negative of J-"^^^^{p) in Eq. ()A2|1 with 6m^ given by Eq. ()A4jl with 
q and p replaced by q and p, respectively. Clearly, one can freely rename the integration 
variable q as q. Also, it is well known that (5m^ does not explicitly depend on p, but only on 

— p^ = m^. Since we have Pq — p'^ = rri^, we may replace p by p. Hence we can conclude 
that J^i^^<^''^s(p) = -J^°^'^'^'^(p). 

Next we analyze the contribution jF<(p). By changing the integration variable from q to 
K = p — q and then to k = K/^, we find 



•^<(P) 



2hpo 



dtidt2 



2go(27r)3 2fc 



X 



^*Jh)Mt2) Vi {UM. P, q)0q(ti)C(t2) 



-ik{t-i—t2) 



(A34) 



Now we have q = p — /ik. (Note also that the upper limit h°'^^\ of integration for k becomes 
infinite as h tends to zero.) Hence, we have q — > p for all k as /i — > because h ■ h"~^X — > 0. 
The exponential factor takes the form 



exp < t 



dp 



2 dp^dp> 



dC \ . (A35) 



Thus, to truncate the series in the exponent at the second term for all k in the integration 
range, we need h{h'^~^)'^ — > as /i — > 0. This is satisfied due to the requirement a > ^ 
Thus we have 



4- 



•^<(P) 



X exp 



dtidt2 
ik{t2 - ti 



d^k 



*^ P - V(C) 



kdc 



[p - V(ti)] ■ [p - Vih)] 



0"p(il)crp(t2) 



^p(C) 



4a— 4^ 



(A36) 



Since 4a— 4 > — 1 because of the requirement a > |, the non-leading terms do not contribute 
to the position shift in the limit h ^ 0. We drop this contribution from now on for this 
reason. Recalling that [p — V(t)]/crp(t) is the velocity of the corresponding classical particle, 
dx./dt, we obtain at leading order, in analogy with Eq. 



d^k 



+ 00 



di 



^ ^' d^ d^' 



(A37) 



where we have defined = ti — n ■ x(ti) and ^' = t2 — n. ■ x(t2) with r\. = k/k. If we write 
the Heaviside function as d{i - i') = 1/2 + - ^')/2, where e(C' - ^) = 1 if C > ^ and 
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e(^' — ^) = —1 if ^' < ^, then the first and second terms give the imaginary and real parts 
of JF<(p), respectively. Thus, twice the imaginary part is 



2ImJ^<(p) 



+00 



+00 



Re^<(p) 



2h 



kdkdQ 
2(27r)3 



(A38) 



This coincides with the emission probability found in Ref. 2| after we let h" — » 00 as 
required by unitarity. The real part is given by 



+00 < n 1 

di'<i'-0-^x{0^x{iY'^''-'^ , (A39) 
00 '^s '^s 



We have introduced the cut-off function x(^) as in Eq. (jl^ . Now, we integrate by parts 
with respect to the variable ^ and add the result of integrating by parts with respect to 
then divide by two. Thus we find 



Re^<(p) = 



dkdVL 



r+oo r+oo f i u 



d d 



di' di) dC 



d^ 



di 



(A40) 



The factor multiplying e*'^'-^'"^^ in the second term is symmetric in ^ and Hence the 
A;-integral of e*'^^^ can be made into S{^' — ^) hy extending the integration range of k to 
(—00, +00) and dividing by two in the limit h 0. The factor inside the square brackets 
becomes zero if one lets C,' = C,- Hence the contribution to Re J-'<(p) from the second term is 
of order higher than h~^. Evaluating the first term with the use of {dx^/ dt) {dx^/ dt) = 1 — 
and ^ = 1 — n • V, we find 



ReJ^<(p) 



2-a 



f 1 - v2 

dt I dVL 

1 — n ■ V 



(A41) 



to order h ^. Now, let us consider jF<'°(p). By substituting the small- approximation 
flA26|) in Eq. ()A31|) and noting the equations |0p(t)P = PQ/a-p{t), <7p(t) = — v^, 

[p - V(t)]/cTp(t) = V and [ap{t)f - [p - V(t)] ^ = we find that J^<'°(p) is indeed equal 
to the leading term of ReJF<(p) given by Eq. ()A41|) . 

Thus, we have shown the equality J^^^^^[p) + jF^'^'^'^™s(p) = and demonstrated that 
jF<(p) — J-'^'^{p) is of order h^^ but is purely imaginary at this order. Therefore, from 
Eq. (jA32j) we conclude that the real part of the forward-scattering amplitude, ReJF(p), is of 
order higher than h~^, and hence does not contribute to the change of position of the scalar 
particle given by Eq. (jTHl). 
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